home *** CD-ROM | disk | FTP | other *** search
/ SGI Hot Mix 17 / Hot Mix 17.iso / HM17_SGI / research / lib / obsolete / simpregress.pro < prev    next >
Text File  |  1997-07-08  |  907b  |  45 lines
  1. ; $Id: simpregress.pro,v 1.2 1997/01/15 04:02:19 ali Exp $
  2. ;
  3. ;  Copyright (c) 1991-1997, Research Systems Inc.  All rights
  4. ;  reserved. Unauthorized reproduction prohibited.
  5.  
  6.  
  7. function simpregress,X,Y,W,YFit,A0,sigma,FTest,R,RMul,Chisqr
  8. Sx = N_Elements(X)
  9.  
  10. if N_Elements(W) EQ 0 THEN BEGIN
  11.   W = Replicate(1.0,Sx)     
  12.   NoWeight=1
  13. ENDIF ELSE NoWeight = 0
  14.  
  15.  
  16.  
  17. Sum = Total(W)
  18. Sumx =Total(W*X)
  19. Sumy = Total(W*Y)
  20. SumX2 = Total(W*X*X)
  21. SumXY = Total(W*X*Y)
  22. SumY2 = Total(W*Y*Y)
  23.  
  24. Del = Sum*SumX2 - Sumx^2
  25. A0 = (sumx2*sumy - sumx*sumxy)/del
  26. B =  ( sumxy * sum - sumx*sumy)/del
  27.  
  28. if NoWeight EQ 1 THEN   $
  29.    Var = (sumy2 + A0*A0 *Sum + B*B*SumX2 - 2.*(A0*sumy + B*sumxy - A0*B*sumx))  $
  30.          /(Sx-2)                $
  31. ELSE Var = 1.0
  32.  
  33.  
  34. sigma = sqrt(Var*Sum/Del)
  35. RMul = (sum * sumxy - sumx * sumy)/sqrt(Del*(sum*sumy2 - sumY^2))
  36. R = RMul
  37.  
  38. YFit = b * X + a0
  39. Resid = Y - YFit
  40. chisqr = Total(W*Resid^2)
  41.  
  42. return,b
  43. END
  44.  
  45.